Stability and pairing in quasi-one-dimensional Bose-Fermi mixtures 
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We consider a mixture of single-component bosonic and fermionic atoms in an array of coupled 
one-dimensional "tubes" . For an attractive Bose-Fermi interaction, we show that the system exhibits 
phase separation instead of the usual collapse. Moreover, above a critical inter-tube hopping, all first- 
order instabilities disappear in both attractive and repulsive mixtures. The possibility of suppressing 
instabilities in this system suggests a route towards the realization of paired phases, including a 
superfluid of p-wave pairs unique to the coupled-tube system, and quantum critical phenomena. 
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Recently, heteronuclear resonances in mixtures of 
bosonic and fermionic ultracold atoms have attracted no- 
ticeable theoretical and experimental interest, due to the 
possibility of generating and exploring novel quantum 
phenomena in a controllable manner. For example, by 
varying the interaction in a Bose-Fermi (BF) mixture, 
one can, in principle, observe a quantum phase transi- 
tion from a Bose-Einstein condensate (BEC) to a nor- 
mal Fermi gas phase by binding bosons and fermions 
into fermionic molecules [l], Q. Indeed, this feature 
has already been exploited to create deeply-bound, polar 
fermionic molecules Q. However, the single biggest im- 
pediment to realizing such novel phenomena in BF mix- 
tures is substantial inelastic collisions. The situation is 
particularly severe on the attractive side of the heteronu- 
clear resonance, where a collapse of the cloud has been 
observed 0, [B| , resulting in a sudden loss of atoms from 
three-body recombination. On the repulsive side of the 
resonance, an interaction-induced spatial separation of 
bosons and fermions 0,0] ensures that the atomic system 
is relatively stable [g, [9(. However, if one sweeps through 
the resonance, the system once again suffers significant 
inelastic losses when molecules collide with atoms fioj ]. 

In this Letter, we argue that many of these obstacles 
may be circumvented by embedding the mixture in a 
two-dimensional (2D) array of ID tubes generated via 
an anisotropic optical lattice. Such a lattice is experi- 
mentally realizable and has alread y b een used to explore 
the 1D-3D crossover in a Bose gas . While strictly ID 
BF mixtures have been investigated extensively in sev 
eral theoretical works [l2, 13|, 14 



we demonstrate using mean-field theory that, similarly 
to ID [l6| and contrary to expectation 00, there is 
no collapse in a quasi-lD attractive mixture. Crucially, 
we find that the hopping can be used to suppress first- 
order instabilities in BF mixtures and, as such, it may al- 
low one to investigate quantum phase transitions induced 
by BF pairing [HQ, without the intrusion of first-order 
transitions. In addition, we will show using the Lut- 
tinger liquid formalism that, for a sufficiently strong BF 
attraction, the coupled-tube system exhibits an exotic 
superfluid phase, where p-wave pairing occurs between 
fermionic molecules comprised of a single boson and a 
single fermion. 

In the following, we consider a mixture of bosonic (b) 
and fermionic (/) atoms confined in an N x x N y square 
array of ID tubes of length L z . We focus on the homo- 
geneous case, but our results can easily be mapped to 
the case of a harmonic trapping potential using the local 
density approximation Q . For sufficiently strong lattice 
confinement, the xy motion can be approximated by a 
single-band, tight-binding model (setting h = 1), 
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+ 2t [2 - cos(k x d) - cos(kyd)} , (1) 



where t is the hopping between tubes and d is the tube 
spacing. Here, the transverse xy momenta are restricted 
to the first Brillouin zone, \k X:V \ < ir/d. The single- 
channel Hamiltonian is thus 



the novelty of our approach is to allow a finite hopping 
between tubes, thus preserving the true long-range order 
of condensed phases as found in 3D, while still main- 
taining the advantages of a ID system. In particular, 3- 
body recombination should be greatly reduced, perhaps 
even more than in a BF mixture confined to a 3D optical 
lattice (see, e.g., (HI), since its rate vanishes for short- 
ranged interactions in the ID limit [22| . Furthermore, 
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where = e(^ b — /Xf^ and //f (//{,) is the fermionic 
(bosonic) chemical potential. The contact interactions 
Ubf, Ubb are effectively ID, and we choose a repulsive 
boson-boson interaction Ubb > to ensure the stability 
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of the Bose gas. Interactions between identical fermions 
can be neglected due to the Pauli exclusion principle. If 
all atoms experience the same transverse trapping fre- 
quency lu±, then the ID interactions, Ubf and Ubb can 
be written simply in terms of the 3D scattering lengths 
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where the oscillator length of the tube a a/ 3± = 
yJl/m a pujL depends on the masses uibb = rrib and 
rriBF = 2mfmb/(mf + mb), while C ~ 1.4603/V2- 

The introduction of an inter-tube hopping t naturally 
leads to a crossover from ID to 3D behavior. The limit 
( \t ^ * recovers the isotropic 3D dispersion, while the 
opposite limit e£' b ^ 8t corresponds to the ID regime. 
For degenerate fermions, this implies 3D behavior when 
the Fermi energy Ep <C t, i.e. at sufficiently small den- 
sities, and ID behavior when ep ^ 8t, i.e. at large 
densities. However, for weakly-interacting degenerate 
bosons, the spread of the momentum distribution is set 
by the temperature T and thus we require ksT <C t and 
ksT ^> 8t, respectively, to access the 3D and ID regimes. 
A corollary of this is that we expect the superfluid criti- 
cal temperature T c of the quasi-lD Bose gas to be finite 
and scale as some positive power of t. Contrast this with 
the strictly ID limit (t — 0), where T c is strictly zero. 
We shall focus on the T = limit, so the effective dimen- 
sionality will only depend on the fermion density. 

We begin by analyzing the first-order instabilities of 
the quasi-lD mixture using mean-field theory. Of course, 
for purely ID mixtures, a mean-field description [l3|, [l6| 
is unreliable because the physics is dominated by fluctu- 
ations, and one must instead use the Luttinger liquid for- 
malism (3]. However, we expect a mean- field treatment 
to be reasonable for finite inter-tube hopping, because 
then it works well in the low-density 3D limit (t/ep ^ 1), 
as well as being consistent with the Luttinger liquid de- 
scription in the high-density, weak-coupling, ID regime 
(8t/e F < 1, \U B f\ y / 2m f /e F < 1) Q- Specifically, we 
take 6k = o~k,o\/ L z N x N y <fr, so that the grand-canonical 
free energy density Q([if, Mb) = min$ /($, /!/, lib) can be 
easily evaluated by integrating out the fermionic degrees 
of freedom, giving: 
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In addition, the ID densities of fermions and bosons 
in each tube are given respectively by rib = & 2 and 
rif = ^/(N x Ny)^2^' z k kpz/ir, so that, within mean- 
field, we always have a BEC when rib > 0. Here, the 
system dimensionality is set by the parameter ([if — 
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FIG. 1: (Color online) Zero temperature, mean-field phase 
diagrams in density space for a repulsive BF mixture with di- 
mensionless hopping strength t' = 2m/t/c| = 0.5 (left panel), 
and an attractive mixture with t' — 0.2 (right). Each mixture 
can either form a uniform BEC phase (light gray shaded re- 
gion) or it can undergo a first-order transition (thick red lines) 
to a phase-separated state (dark gray). The dotted lines con- 
nect points on the first-order boundary with the same chem- 
ical potential. Filled circles mark stable tricritical points, 
where 1 st and 2 nd (rib = 0) order transition lines merge, while 
the empty circles correspond to unstable ones. Spinodal lines 
(blue dashed) divide the phase-separated region into the un- 
stable domain (internal region) and the metastable domain 
(external). Note that, for the repulsive mixture, the phase at 
very low fermionic densities is always uniform. 



UBF& 2 )/t or, equivalently, irrif / a^2trrif , where a = 
Jq dk x dk y ^2 + cos k x + cos k y /ir 2 ~ 1.35. 

For a 3D, attractive (Ubf < 0) mixture with no optical 
lattice, it is easy to see that the free energy at large $ 
is dominated by the BF interaction term (cx — $ 5 in 3D) 
and is thus not bounded from below [27| . This implies 
that the system is unstable to collapse at sufficiently high 
densities [7|, |24|. On the other hand, in a ID tube, the 
interaction term instead scales like — $ 3 at large <& and 
is thus compensated by the boson-boson repulsion (cx 
4> 4 ) [l6| . Therefore, contrary to what has been previously 
assumed both ID and quasi-lD attractive BF 

mixtures will exhibit phase separation instead of collapse. 

After minimizing the free energy ([4]) with respect 
to the boson field <f>, we can construct the phase di- 
agram using just three dimensionless parameters, such 
(4) as the dimensionless hopping strength t' = 2mft/c 2 , 
and the dimensionless densities ribj/cbj, where q, = 
2m f \U B F\ 3 /(K 2 U BB ), c f = 2m f U B p/(ir 2 U B B)- The re- 
pulsive, strictly ID (t = 0) case has been evaluated 
within mean- field in Ref. Here, contrary to the 3D 
case 0], phase separation occurs at low fermionic den- 
sities, rif/cf < 3/4, irrespective of the boson density. 
Furthermore, phase separation only occurs between two 
pure phases (when rib/cb < 3/4) or between a mixed 
phase and a purely bosonic phase (when rib/c-b > 3/4). 
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FIG. 2: (Color online) Evolution of the stable (filled blue cir- 
cles) and unstable (empty circles) tricritical points n^ cp /c/ 
( n & CP = 0) as a function of the hopping strength t' for a re- 
pulsive mixture. For attractive mixtures, the tricritical points 
have the same values, but the stable and unstable branches 
are switched. The (dashed) line approximately separating the 
ID from the 3D regime is rif/cf = \/2Fa/iv, with a ~ 1.35 
— this defines the point, /if — Ubf& 2 = 8t, where the Fermi 
surface first touches the xy band edge. Tricritical points dis- 
appear altogether above the critical value t' CT ~ 1.19. 



The topology of the repulsive phase diagram changes sub- 
stantially once t > 0. As shown in Fig. [TJ a stable tri- 
critical point [28| appears at low fermionic densities and 
there is instead a uniform phase for n/ / \J1tmf <C a/ it. 
Here, the phase diagram resembles the 3D phase diagram 
derived in Ref. 0], as expected. By contrast, at higher 
fermionic densities, we recover ID behavior, such that 
phase separation only exists for nj/cf < const. However, 
we note that there is never phase separation between two 
pure phases at finite t, unlike in the strictly ID and 3D 
cases. Instead, phase separation either occurs between 
a purely fermionic and a mixed phase (for nj/cf < 0.73 
and rib/cb < 0.73 in Fig. |T|) or between two mixed phases. 

For the attractive case, the structure of the phase di- 
agram does not change when hopping is switched on. 
Moreover, unlike the 3D case, the mixture displays phase 
separation instead of collapse, as previously discussed. 
However, the phase diagram has a region where phase 
separation occurs between a mixed phase and the vac- 
uum (see Fig. [T]), and this may be viewed as a remnant 
of the collapse in the 3D system. Note that the tricriti- 
cal points have the same values as in a repulsive mixture 
at the same t', but their stability is switched. Both at- 
tractive and repulsive mixtures at the same t' also feature 
identical spinodal lines [29I ] , which indicate when the sys- 
tem becomes linearly unstable to phase separation. 

By tracking the evolution of the stable and unstable 
tricritical points as a function of t' (Fig. [2|), we find that 
the situation dramatically changes at larger t' . Notably, 
at the critical value t' CI ~ 1.19, the stable and unstable 
tricritical points merge, and the width of the phase sep- 



arated region reduces to zero. Thus, for t' > t' CT , the sys- 
tem exhibits only a uniform phase. This is a consequence 
of the 1D-3D crossover in this system: eventually the in- 
stabilities of the ID regime fall in the low density regime 
where 3D behavior dominates, and vice versa. These re- 
sults suggest that one can stabilize a BF mixture using an 
appropriate 2D optical lattice. Indeed, we find that if the 
oscillator length o,bb± is comparable to the boson-boson 
scattering cibb, then we can have the situation where 
there is phase separation or collapse in 3D and yet no 
instabilities in the quasi-lD case for t > t' cr . This is ba- 
sically because the effective ID interaction Ubb diverges 
at asB± = Cclbb in Eq. j3]). 

The absence of collapse makes quasi- ID systems ideal 
for examining BF pairing. In particular, the possibility of 
suppressing first-order instabilities by tuning the hopping 
strength opens up the prospect of investigating quantum 
phase transitions. For example, one could realize a con- 
tinuous quantum phase transition, where the BEC is de- 
stroyed by the pairing. Even though our mean-field treat- 
ment does not include the possibility of pairing, we can 
still get an estimate of the instability towards BF pairing 
by considering the two-body problem. In the t — limit, 
the BF binding energy is approximately e = ?71bf£/r F /4, 
which becomes exact in the limit \clbf\ — * [25(. How- 
ever, once we switch on the hopping, the bound state is 
lost when t ~ 0.05rriBFU BF . Thus, if we wish to explore 
pairing-induced phase transitions in the absence of first- 
order instabilities, we need this "resonance" to lie above 
t' cr , i.e. we require Ubb/Ubf ^ 0.2(1 + mf/rrib). 

We can determine what symmetry-broken states may 
exist in the quasi-lD system by comparing the decay of 
different correlation functions in the purely ID limit. On 
general grounds we expect that the operator with the 
slowest decay in ID will fix the long-range ordering of 
the higher-dimensional system. This is due to the fact 
that the exponent 77 governing the spatial decay of the 
operator O, (O x Oo) ~ also appears in the suscep- 

tibility as a function of temperature, x(T) ~ T r >~ 2 , and 
one can show using a mean-field approximation for the 
intertube couplings that the symmetry-broken state with 
the most divergent x(T) generally has the highest T c [2^ . 

In the strictly ID limit, the low-energy low-wavelength 
effective field theory is described by the Luttingcr formal- 
ism (bosonization) In particular, we consider the 
case where the fermionic and bosonic phase velocities are 
similar Vf = Vb = v and the low-energy effective Hamil- 
tonian can be described by introducing in- and out-of- 
phase phase and density fluctuations of the mixture 1J] , 
0i,2 = 4^(06 ± 0/) and 0i, 2 = -j=(6 b ± 9 f ) : 
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where AA2 are Luttingcr parameters and A a cut-off. 
The in-phase mode 1 describes a one-component (gap- 
less) Luttinger liquid. We emphasize that Eq. ([5]) holds 
for any value of the BF coupling. However, ATi,2 can only 
be determined analytically in the limit of small Ubf and 
one must resort to numerics [20[ when a perturbative ex- 
pansion in Ubf is no longer accurate. 

In the limit of equal filling, 71/ = rib, the field $2 ac- 
quires a gap — the corresponding "paired" phase has 
been introduced in Ref. 14J . The slowest algebraic decay 
is then given by the operator O — bf when K\ > 2/V3; 
otherwise it's given by charge-density wave correlations. 
However O is a fermionic operator and as such cannot 
lead to condensation when we couple the tubes to access 
the 3D limit. Instead we must consider the composite 
bosonic operator C>(") — whose correlations can 

be evaluated from Eq. (0: 
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In the paired phase the field is pinned by the relevance 
of the cosine operator in Eq. ([5]) and therefore the con- 
jugate field 62 has exponentially decaying correlations. 
As a consequence, <D^ n ' also has exponentially decaying 

correlations unless n = 2. Eq. ((6|) leads to a decay law 

~ (2) A (2) 4 
(Ox Oq ) ~ l/x K i and dominates over charge-density 

wave correlations for K\ > 2. Therefore, in a system of 
weakly coupled tubes we expect condensation of the opera- 
tor <D^ 2 \ i.e. a p- wave paired phase of fermionic molecules 
each comprised of a single boson and a single fermion. 
In the strictly ID system, this operator was recognised 
to have quasi-long-range order in Ref. [3]. Here, we 
find that, for K\ > 2, the operator O^ 2 ' has the slow- 
est decaying correlations and therefore implies condensa- 
tion of fermionic molecules in the higher-dimensional sys- 
tem. Moreover, the p-wave phase of fermionic molecules 
is topologically distinct from a superfluid of p-wave pairs 
of atomic fermions coexisting with a BEC. In the for- 
mer, even though the global phase symmetry is bro- 
ken, there is a remaining subgroup U(1)b-f of rela- 
tive phase transformations between bosons and fermions 
that is preserved, since b 2 and JlJr can be rotated by 
opposite phase factors without changing the superfluid 
order parameter. The p-wave order parameter Ak also 
breaks full spatial rotation symmetry 50(3) and, in the 
ground state, can either be of the form Ak <x k z (a spin- 
less equivalent of the polar phase of superfluid 3 He) or 
Ak oc k x + ik y (the spinless variant of the 3 He A-phase). 
Finally, we note that it should be possible to access this 
p-wave phase in an attractive BF mixture: a numerical 
DMRG analysis for 40 K- 87 Rb mixtures [U has shown 
that deep inside the paired phase there is at least one 



point with K\ ~ 2. 
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<9^2)2/|$ = 0- 



/Xi,$ 2 + Ubb$ 4 /2, with k F = k Fz (t 
9 { 3 s2)3 /|*=o > and d&i /|*=o = 
At the spinodal lines 8^,2/ = = d^2j2f = 0. For a 
given t' the fermionic density along each spinodal line is 
fixed at the density of a tricritical point, n/ = n^ p . 



